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INTRODUCTION 


Many future NASA missions would utilize significantly large and 
flexible spacecrafts and would require very stringent pointing and 
vibration suppression requirements. The active controller that can 
achieve these objectives will have to be designed with very accurate 
knowledge of the dynamic behavior of the spacecraft to ensure per- 
formance robustness to a variety of disturbances and uncertainties. 

In the past few years, several design approaches were proposed 
for vibration control during and after slewing maneuvers. NASA 
Langley Research Center initiated the Spacecraft Control Laboratory 
Experiment (SCOLE) program [1] to promote direct comparison and 
realistic test of various control design techniques against a common 
laboratory article. The article was intended to resemble a large 
space antenna attached to the space shuttle orbiter by a long 
flexible mast. 

The primary control objective of SCOLE is to direct the RF line- 
of-sight (LOS) of the antenna-like configuration towards a fixed 
target under conditions of minimum time and limited control 
authority. 

This problem of directing the LOS of antenna -like configuration 
is studied as being composed of two-phase control problem during the 
current research period. The two phases are namely, the slew 
maneuver control of rigidized-body of SCOLE configuration and the 
vibration suppression of flexible antenna. This formulation allows 
the design of control systems using a decentralized control scheme 
in which the dynamics of the two phases of control problem are 
viewed as two subsystems with some interaction. The residual 
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vibration suppression of the flexible antenna at the end of the slew 
maneuver is viewed as the second phase of the control problem. 


BRIEF SUMMARY OF MAJOR ACCOMPLISHMENTS 

(a) Decentralized Slew Maneuver of SCOLE 

During the current research period, the basic software for the 
decentralized slew maneuver control of SCOLE model was completed. 

The software was based on an algorithm of Hierarchical Optimal 
Control formulation in which the SCOLE model is viewed to be a 
large-scale dynamical system composed of interconnected subsystems. 
It was deduced that the system could naturally be decomposed along 
the lines of dynamics of rigid part and those of flexible antenna. 
Thus the set of uncoupled subsystems was generated together with 
their coupling relations. The performance index to be minimized was 
written in terms of the two subsystems and the subsystems were 
further formulated individually in terms of two-point boundary value 
problems. The necessary conditions for the individual subsystems 
were developed in terms of nonlinear differential eguations for 


optimal control schemes. 

(b) Quasilinearization of Subsystems: 

The two-point boundary value problem of each subsystem was solved 
by the method of quasilinearization. The process of quasilineariza- 
tion for each subsystem involved the linearization of the nonlinear 
system equations and utilizing the method of complementary 
functions. An iterative search was incorporated in the algorithm to 
get the final set of equations which satisfy the split boundary 
conditions . 
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The complete control system design was developed at a second 
level in terms of a coordinating algorithm to optimize the overall 
state trajectory. 

The development of the software to design an optimal state 
feedback control system for each subsystem allows specifications of 
arbitrary large-angle nonlinear slew maneuvers. Computer simula- 
tions of second-level trajectory optimizations were obtained and are 
being analyzed for determining global optimal solutions. 

(c) Shooting Method 

The two-point boundary value problem for each subsystem is also 
being solved using shooting method instead of quasilinearization 
method to compare the solutions in terms of numerical convergence. 

(d) Presentations and Publications 

A detailed technical report based on Slew Maneuver Control has 
been prepared for publication as NASA Contractor Report. A 
presentation entitled, "Decentralized Slew Maneuver Control and 
Vibration Suppression of Large Flexible Spacecrafts," was presented 
at the Workshop on Computational Techniques in Identification and 
Control of Flexible Flight Structures, Lake Arrowhead, California, 
November, 1989. This paper is also to be published in the 
proceedings. A second paper entitled, "Slew Maneuvers of Large 
Flexible Spacecrafts" is to be presented and published at the 1990 
American Control Conference in May 1989. 

A chapter based on collective work of previous years is to be 
published in a book entitled, "Modelling and Control of Large Space 
Structures: The SCOLE Experience" to be edited by A. V. Balakrishnan 
and L. Taylor. This is included in the appendix. 
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SLEW MANEUVER DYNAMICS AND CONTROL OF 
SPACECRAFT CONTROL LABORATORY EXPERIMENT 

(SCOLE) 


Y. P. Kakad 

Dept, of Electrical Engineering 
University of North Carolina at Charlotte 
Charlotte, NC 28223 


ABSTRACT 

In this article, the dynamics and control of slewing maneuvers of NASA 
Spacecraft COntrol Laboratory Experiment (SCOLE) are analyzed. The control 
problem of slewing maneuvers of SCOLE is formulated in terms of an arbitrary 
maneuver about any given axis. The control system is developed for the combined 
problem of rigid-body slew maneuver and vibration suppression of the flexible 
appendage. The control problem formulation incorporates the nonlinear dynamical 
equations derived previously in a report [l] and is expressed in terms of a two- 
point boundary value problem utilizing a quadratic type of performance index. 

The two-point boundary value problem is solved as a hierarchical control 
problem with the overall system being split in terms of two subsystems, namely 
the slewing of the entire assembly and the vibration suppression of the flexible 
antenna. The coupling variables between the two dynamical subsystems are 
identified and these two subsystems for control purposes are treated independently 
in parallel at the first level. Then the state-space trajectory of the combined prob- 
lem is optimized at the second level. 


This work was supported by NASA grant NAG-1-535. 
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1. INTRODUCTION 

The primary control objective of the Spacecraft Control Laboratory Experi- 
ment (SCOLE) is to direct the RF Line-Of -Sight (LOS) of the antenna-like 
configuration towards a fixed target under the conditions of minimum time and 
limited control authority [2]. This problem of directing the LOS of antenna- like 
configuration involves both the slewing maneuver of the entire assembly and the 
vibration suppression of the flexible antenna-like beam. The study of ordinary 
rigid-body slew maneuvers has received considerable attention in the literature 
[3,4] due to the fact that any arbitrary large-angle slew maneuver involves 
kinematic nonlinearities. This is further complicated in the case of SCOLE by vir- 
tue of a flexible appendage deployed from the rigid space shuttle. The dynamics of 
arbitrary large-angle slew maneuvers of SCOLE model were derived previously by 
the author in a report [l 3 as a set of coupled equations with the rigid-body motions 
including the nonlinear kinematics and the vibratory equations of the flexible 
appendage. These nonlinear and coupled dynamical equations are used in this arti- 
cle to study the slew maneuver control in terms of a hierarchical feedback control 
scheme. 

The control problem of slewing maneuvers of this large flexible spacecraft is 
developed by using the two-point boundary value problem in terms of the rigid- 
body slewing and the vibration suppression of the flexible appendage as two 
separate dynamical subsystems. A decentralized optimal control scheme is utilized 
in order to solve individual boundary-value problem for each of the two subsys- 
tems by defining their state variable models and incorporating the coupling vari- 
ables between the two subsystems in these models. Also, the boundary conditions 
of the overall system are reworked in terms of boundary conditions of each sub- 
system. A quadratic performance index is utilized for the overall system and is 
subsequently expressed in terms of a sum of two individual performance indices of 
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the subsystems. 

The basic algorithm for obtaining an optimal closed-loop state feedback 
scheme involves using a trajectory in terms of a vector of Lagrange multipliers as 
an initial guess at level two. This is used at level one in quasilinearization applica- 
tion. 

The two-point boundary value problem for each subsystem is solved at level 
one by using a quasi linearization technique as a trajectory optimization problem. In 
the quasilinearization procedure, starting from an initial guessed state trajectory, 
successive linearizations are performed of state equations in such a way that the 
final solution of the state trajectory is within an acceptable degree subject to the 
boundary conditions. The state vector definition at this level is an augmented state 
vector which includes both system states and costates. 

These optimum solutions of the subsystem trajectories are utilized at level 
two to yield the updated trajectory of the vector of Lagrange multipliers of the 
overall system to be used for quasilinearization process at level one. The basic 
steps of the algorithm are repeated to optimize this second level trajectory with 
respect to a prespecified error criterion to obtain an optimal feedback law. 


-4 - 


a_ 

B 

F 2 {t) 
Go(? ) 
h 

li 

J 

K 

L 

M 

N 

a 

<h 


r 


X 


r 


y 


Mi 

3k 

Zi 

X 


$ 'xl 
<l>yi 


2. LIST OF SYMBOLS 

Vector to the point of force application on the beam 
Damping matrix 

Force applied at the reflector mass center 
Moment applied about the orbiter mass center 
Equivalent Mass moment of inertia of total assembly 
Mass moment of inertia matrix of the reflector 
Functional used for two-point boundary value problem 
The stiff ness matrix 
The Length of the beam 
Angular velocity vector transformation 
Effective moment applied at the reflector c.g. 

The total number of subsystems 
The generalized force vector 
Generalized coordinates 

Position vector f rom the orbiter mass center to the point of 
attachment 

x co-ordinate of the reflector mass center in the body-fixed 
frame 

y co-ordinate of the reflector mass center in the body-fixed 
frame 

Control vector of i th system 

State vector of i th system 

Vector of interconnecting variables 

Unit vector representing the axis of rotation during the 
slew maneuver 

i th Eigenfunction corresponding to u x 
i th Eigenfunction corresponding to Uy 
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i th Eigenfunction corresponding to 
The attitude of the orbiter in the inertial frame 
Slew Angle 

The angular velocity of the orbiter in the inertial 
frame 

The angular velocity of the reflector in the inertial 
frame 

Vector of Euler parameters 
j ) Direc delta f unction 

Dual functional for two-point boundary value problem 
Vector of Lagrange multipliers 
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3. ANALYTICS 

Slew Maneuver Specification and Control Variables 

The analytics for the dynamics of SCOLE developed in reference [l] are used 
to derive the control laws for an arbitrary slew maneuver. It is assumed that the 
slew maneuver is performed by applying moments on the rigid shuttle and the 
vibration suppression is achieved by means of forces on the flexible antenna and 
the reflector. The slew maneuver is considered to be an arbitrary maneuver about 
any given axis [l]. 

The slew maneuver is defined in terms of x the axis about which the slew 
maneuver is performed. If £ is the slew angle and <o is the angular velocity of the 
orbiter in the inertial frame, then the four Euler parameters can be defined as 


€1 = yisin-^- 

e 2 = y 2 sin-| 

= y 3 sin-| 

£ 

€4 = COS-f- 
2 


Cl) 


The four Euler parameters can be related to the angular velocity components 
of the rigid assembly as 


*1 


€ 1 *4 — *3 *2 


0 

«2 


e 2 *3 U ~ C 1 


6>1 



€3 “€ 2 «1 *4 


0>2 

*4 


€4 €j — e 2 -e 3 


«> 3 


The slewing maneuver can be given in terms of the following equations [l] 


I 0 £> + A £ = Q_{t )+££ 2 (<o) 


(3) 
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A 2 T u + Ail + 5i. + K± = Q{t) (4) 

where, 

0.(7 ) is the net moment applied about the mass center of the orbiter and is 
given by the following equation (figs. 1 & 2) 

a(r) = &(0 + (r +o.)x£ 2 . (5) 

Also, ) represents the generalized force vector which is given by the following 
equation 


acr) = 


Z ( Qjxfc ) + Qjyfj )) + Qx, + Qyj + Qv», 

j =i 

Z C Qj Xa it ) + Qjyp )) + Q Xi + Q yi + Q + 

)• i 


Z c Qjx t b ) + QjyS* ^ + Qxi + Qyi + Q*i 

7-1 


( 6 ) 


where, the generalized f orce components are given as 


L 


Qjx t = J F JX C Z ~ Z j toxi tez 

0 
L 

(7) 

Qjy t = f Fjy(z ,t )8(z- 2; )</>yi(r )dz 
0 

OO 
' ' 

O 

/ — \ 
* 

II 

o 

• 

(9) 


Here, F jx (z,t ) is the x component of the concentrated force applied at location j 
on the flexible antenna and Fjy is the y component of that force. 
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Also, 


Q*^) = F 2x (r )0 xi U) 

Qyi(t) = F 2 y tT)<t> yi U) ( 10 ) 

<2^(0 = M v& (r)0^(L) . 

Here, F 2 is the force applied at the reflector C. G. 

Thus, 

) = F^Ty + F ly r x + M 2 ^ (11) 

The location of reflector C. G. is given by coordinates (r x ,r y ) and M 2x p 
represents the external moment applied at the reflector C. G. Also, the nonlineari- 
ties N 2 can be expressed in terms of pure rigid body kinematic nonlinearity and 
the nonlinear coupling term between the rigid-body modes and the flexible modes. 

N 2 = A 4 ( co,9_ ) + A 5 ( (12) 

The details of this term N 2 are included in appendix A and the equation (12) 
can be further simplified in terms of Euler parameters by relationships developed 
in Appendix B as 


N 2 = A 6 (o>,_i) + A 7 ( <&€_)£. (13) 

where e_ is the Euler vector comprising all four Euler parameters. 

From equations (3) and (4) and by defining A = — A 2 7 0 -1 A 2 + A 3 , the 
following equations are obtained 


<£. = I o 1 


A 2 A ~ l Bg_ + A 2 A ~ l Kg_ + 


U 2 A - 1 A 2 r / e - 1 +/ 3 


a(r) 


+ U 2 A- l A 2 r I 0 - 1 +I 3 \ 




(14) 
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l = -A-'Bi - A - 1 Kg. - A- x Afl- x QXt ) + A~ l A 2 T 2 (tt**) 

+ A~ x a Cf ) • (15) 

It is assumed that control forces applied for vibration suppression has negligi- 
ble effect on rotational maneuver of the spacecraft in developing equations (14) 
and (15). Also, / 3 represents 3x3 identity matrix in these equations. 

Subsystems and State Variable Models 

The two dynamical subsystems considered for decentralized control are the 
dynamics of the slewing of the rigidized SCOLE assembly and the vibration 
dynamics of the flexible antenna. These subsystems are represented by subscripts I 
and II respectively for subsequent analysis. 

The following are the definitions of state variables and control variables for 
subsystem I. 


x/iAfi; *i 2 A ei' x i3 A x 1A a c 4 ; 
x I6 Ao 2 ; x n A<* 3 ‘> 

XjsAGu x i9 A G 2» *noAG 3 ; 

u Il AG l ; U; 2 AG 2 ; U/3AG3 . 

The interconnecting variables from the second subsystem to reflect the cou- 
pling between the subsystems are defined as 

z Il A X IIU z I2A x JI 2> z IiA x II3> Zj A A x II4> z IsA®'> 

Z 16 AO ; Z/ 7 AO ; 

Zj 3 A 0 1 Z I9 A ^ I 
2/loAO . 
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The following state equations are obtained for subsystem I using these 
definitions 


. 




*7 1 


*71 


*72 


*72 


*73 


*73 


*74 


Oi D ! 0 


*74 


0 1 0 l 0 



- r ~ f -- 


— 


— 4- — H 

X I5 


1 


X I5 


1 

1 



a: 

O 

o 


X I6 

+ 

0 1 b 2 ' b , 

X 1 6 


1 l_ 



--t- - 

! *• 
1 


1 1 
0, 0 , 0 


x n_ 


O 

o 

*7 8 


X IS 





X I 9 


*79 





*710 


*710 





“/i 
U I2 
“73 
2/1 
2 1 2 
2 12 
2 1 4 


+ 


0 

H 

0 


A^2 ( ^7 1 >*/ 2 >*7 3 >*7 4 >*7 5 >*7 6 >*7 7 ) 


Here.H = / 0 _1 [ A 2 A -1 A 2 r / 0 _1 + / 3 , 


b 2 = i 0 - 1 a 2 a~ 1 k , 


B % = I Q ~ l A 2 A~ l B , 


( 16 ) 


Z) = 


*74 — *73 *72 

*73 *7 4 — *71 

— *72 *71 *74 

“ *71 ~ *7 2 — *73 


For the second subsystem which is the flexible appendage of the entire system, 


the first two flexible modes are considered and the corresponding state variables 
and control variables are defined as 
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x II2&<l2> 

*m Aiu 

u iii AQv umAQti 

As in the previous case, the coupling between the two subsystems is derived 
in terms of the following interconnecting variables from the first subsystem. 


2 u\ A x i\\ 

2 II 2 A X 1 2 > 

2 II 3 A x I3’ 

2 iiaA x ia> 



2 IIS^. X IS> 

Z II6 A x I6> 

2 iii A x n'> 

2 in A*/8l 

2 119 A x I9\ 

Z 11 10 A *710 • 


The following are the state equations of this subsystem. 


x n 1 
X II 2 
x n 3 
X 1I A 



x m 

X II2 
X I1 3 
X IIA 


o I o 


A" 1 | -A~ 1 A 2 t I 0 -' 


u iii 

Ull 2 
2 II 8 
2 II 9 
2 II 10 


0 


-a~ 1 a 2 t i 0 ~ 1 


N 2 ( 2 II 1 > Z II 2< Z II 3 < Z II A < 2 IIS> z II6’ 2 II 7 ) 


(17) 


The Optimal Control Problem 


A general problem for the optimal control of interconnected dynamical sys- 
tems like large flexible spacecrafts can be formulated as 


Minimize J ( Jy & ) i = l,2,...,N 


(18) 
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w.rJ. Uj 

where is the n t dimensional state vector of the ith subsystem, is the 
corresponding m i dimensional control vector and is the r i dimensional vector of 
interconnection inputs from the other subsystem. The integer N represents the 
total number of subsystems and the scalar f unctional J is defined by 


N 

•/ = Z 

i =1 


/ 

Pi ( Cf/ ) ) + f Lt 


i (r 0 ),£ (r ) 


dt 


(19) 


where L t x< (r (r )^ (r ) 


is the performance index at time t for i = 1,2,..,N 

subsystems. The functional J defined in equation (19) Is to be minimized subject to 
the constraints which define the subsystem dynamics, i.e. 


£i = Id 


Xi (r )j^ (r (r ),r 


t 0 <r 7 


x < (t 0 ) = i 0 , i = l,2...,// . 

Also, the minimum of J must satisfy the interconnection relationship 


(20) 


z g/(^ )=a . (2D 

!=1 

The Open-loop Hierarchical Control 

Using the method of Goal Coordination or infeasible method [15,20], we con- 
sider another problem which is obtained by maximizing the dual function <£( X. ) 
with respect to )(f 0 ^ r ^ tj ), where 


<K L . ) = Min 


J (x.JLwZ.A) 


( 22 ) 
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subject to constraints in equations (20) and (21). Here 


*1 


ILi 


1.1 

• 


• 


• 

• 

— 

• 

z_ = 

* 

• 


• 


* 

*/v 


Hn 




( 23 ) 


Also, X^in equation (22) is a vector of Lagrange multipliers which is given as 


ki 


k = 




(24) 


J = Z 

i =1 


Pi ( i (?/ ) ) + f L t ( & Mi xZi ,t ) dt + JVa ( 2L f. ) dt 

t 0 to 

j = 1,2 N 


(25) 


Rewriting this functional J as 


„ N 

J = Z Ji 

i= 1 


N 

= Z 

i = l 


/ 

(i( ?/)) + / Z.i(x i4 Mi l 7i,r) + X r 2i( 

t n 


i ( x, & ,r ) tfr 


(26) 


where, 


V . . 

/i = Pi ( Zi ( tf ) ) + J I L t ( 2 U iMi & J ) + A. r Si ( ^i ^i ,? ) dr (27) 

*0 

and where A. r Gi * ( ^ ,f ) has been refactored into the form X. r Gi ( f ), i.e. 
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into a form separable in the index i . 
Thus, 



( 28 ) 


Then by the fundamental theorem of strong Lagrange duality [22] 


min J = max <J> ( ) , i = 1,2 N. (29) 

Mi k 

Thus an alternative way of optimizing J is to maximize 0 ( X. ) . 

From equation (26), for a given \.( t ) , t 0 ^ t , the functional J is 

separable into N independent minimization problems, the ith of which is given by 


Min 


v - 

Ji — Pi (i ( ) ) + f \ Li ( & ^ & ) + X. r £, (.XixZi ) 

to 1 


dt 


(30) 


X, .Uj , 2 , 


subject to 


ki = h ( 2U Mi £i ) . ?o 
Xi (t 0 ) =x^ 


(31) 


This leads to a two-level optimization structure where on the first level, for 
given k the N independent minimization problems described in equations (30) and 
(31) are solved and on the second level, the X.(i)(r 0 ) trajectory is 

improved by an optimization scheme like the steepest ascent method, i.e. from 
iteration j to j + 1 
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A.(r V +1 = X.(r y + oc j + dj (r) t 0 <r, (32) 

where 

dj = V<KA.CO) = (**) . (33) 

t=i 

V<D( A) is the gradient of $ ( A) . >Ois the step length and dj is the steepest 
ascent search direction. At the optimum dj —* 0 and the appropriate Lagrange mul- 
tipier, X. . is the optimum one. 

The development of this algorithm depends on the assertion Max <£( A) = 
min J and this may not be valid for all nonlinear systems. Consequently, lineari- 
zation of Qj , and linearized equations for Aj may be required for constraints to be 
convex and convexity of the constraints is necessary to prove this assertion. 
Nevertheless, the method is attractive from the standpoint of simplicity and that 
the dual function is concave for this nonlinear case. This ensures that if the duality 
assertion is valid, the optimum obtained is the Global Optimum. 

On the first level, since equation (30) is to be minimized subject to equation 
(31), the necessary conditions lead to a two point boundary value problem from 
which an open loop optimum control could be calculated. However, it is desirable 
to calculate a closed loop control and for this the quasilinearization approach can 
be utilized at level one for all subsystems. Thus an iterative scheme can be set up 
whereby an initial trajectory of A. (O’ , t Q is guessed at level two and 

provided to level one. At level one the two-point boundary value problems of the 
subsystems are solved by quasilinearization. The state and control trajectories of 
all the subsystems obtained at level one are sent to level two. The test for 
optimality based on equation (33) is conducted at level two and if this is not 
satisfied, equation (32) is used to obtain the new A.(f ) for the next iteration. 


Subsystem Closed Loon Controllers 
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The closed loop controllers are obtained at the first level by solving the two- 
point boundary value problems of the subsystems utilizing the quasilinearization 
procedure. As noted in equations (30) and (31), the first level problem for the ith 
subsystem is 

For given A,(r ), t 0 <r , 


min 


Pi 



) + Qi ( 2k ) 


dt 


2k4k>li 


subject to 


(30) 


id = Ld ( 2k 4k «)> U 

2k ( r o ) = 2ko • 

For this problem, the Hamiltonian can be written as 


(31) 


Hi = L t i 2k 4k£i ) + ATQi ) + 1 1i T Ld (2k4k£i ) • (34) 

For a given X, the state and costate equations become 


i» (f ) = Xj (2k 4k & ) 


(35) 


m (r ) = - 


6 Hi 

a 2k 


$k-+&L- k+ 

b2k 8 2 k 



(36) 


with 


a Ui 


= H; 


a^ 

bid 


= 0. 


(37) 


It is assumed here that using the equations (36) and (37), it is possible to 
obtain the control m and the interconnect variable vector Zi which is an explicit 
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function of ru and , i.e. 


Hi = £i ( ,j2i ) 


( 38 ) 


li = dj (**,:& ) 

Using these relationships for u, and z± in equations (35) and (36), the following 
equations are obtained 


Xi = a± ( Xj ) . 

t 0 

(39) 

Hi ^ f ( x tHi ) » 

o 

/A 

/A 

(40) 


with the boundary conditions 


2Li ( to ) 2Li< 

and from the transversality conditions 


(41) 


Hi ( ?/ ) = 


BPi 


xu ) 


Si 


Ouasihnearization Procedure 


(42) 


The two-point boundary value problem of ith subsystem is given by equa- 
tions (39) and (40) subject to boundary conditions of equations (41) and (42). 
This problem is solved by quasilinearization technique as follows. 


Defining y_ = 


Hu 


equations (39) and (40) can be rewritten as 


£(t ) = F_ y_{t ) 


(43) 
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In the quasilinearization procedure, starting from an initial guessed trajectory 
for x = (f ) , successive linearizations are performed of equation (43 ) in such a 
way that the final linear equation for x solves equation (43) to an acceptable 
degree subject to boundary conditions (41) and (42) which could be expressed in a 
more general form as 



X ( t 0 ) T A 0 =b a T 

(44) 


X ( t } y A f =bf T 

(45) 

where A 0 , Af are 2 n x n 

matrices. 


The linearized equation of (43) about a trajectory x = X^ 

(r ) is obtained by 

Taylor series expansion as 



X = 

F (yj) + J (yj)(.y_— yJ)+W 

(46) 


where J {yj ) is the Jacobian of £jx( t ) j, t 0 ^ t ^tf , at yj and 2? represents the 
contribution of the higher order terms. Neglecting these higher order terms, the fol- 
lowing linear equation is obtained 


X = £ ( x' ) + / ( x' ) ( X “ ) • (47 ) 

If the initial guessed trajectory yj while satisfying equations (44), (45) and 
(47) does not satisfy equation (43), then an iterative search can be utilized to 
obtain a better linearizing trajectory by various methods discussed in references 19, 
20, and 21. This iterative search is given by noting that equation (47) can be writ- 
ten by expanding individual equations (39) and (40) by Taylor series expansion 
about a known trajectory xj{t ), r^Jit ),t e , and retaining terms of up to 

first order. The linearized reduced differential equations are 
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i, ^ +1) = a i ^x i j it ),m J it ) 


da 


-ix/ it) ,2/ it)) 


*^ +1) it)-xj it) 


-^- ( x/ it) ,21/ (r ) ) 

cm 


j 

[np +1) Cr ) - 22 / 0 )] (48) 


ib 


O+i) _ 


= bj 


Xi J it ),TU J it ) | + 

8& 


p-ix/it),nj>it)) 

o2L 


xP * 0 it) — xj it): 


a Hi 


(a/CrWCr)) 


2b 0+1) 


(49) 


These differential equations can be rewritten as 


1 0 + 1) = A u it )X, {J+I \t ) + A nit )m 0+1 \t ) + ej it ) (50) 

1 1 0 +1) = A 2l (r )i +1) (r ) + A 22 it ) a ( ' +1) (f ) + e^ (?) (51 ) 
or, in the partitioned matrix form, 


i ( ; + 1 ) (r) 


A ) A 12^ ) 


x/' + 1 ) (r) 


JLi 'CO 

m (j+ 1 \t) 


A 2 i(f ) A 22b ) 


2 h 0 + 1 ) (r) 

+ 

<u J it) 


where the matrices 


(52) 


A nit ) A 


da 


A 12 it ) A 


da 

dUi 


A 2 i(t)A^-, A 22 it) 

fa dm 


and 


£_i J A — A n (r )Xi j it ) — A n it )m J it ) + 


e 2 J A~ A 21 (f )Xi * it ) — A 22 it )m J it ) + A 


( 53 ) 


are evaluated at Xj * it ), 2b J (f ) and hence are known functions of time. 
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The method of complementary functions [24] can be incorporated with this 
linearization of the differential equations in the implementation of iterative search. 

An initial guess, x^° , ,t € [ t 0 , tf ], is used to evaluate matrices in equa- 
tions (53) at the beginning of the first iteration. In the next step, n sets of solu- 
tions to the 2n homogeneous differential equations 


i +1) = A n (t )x t +1) (r ) + A 12 (t ):& +1) (r ) (54) 

ili 0 +1) = A 21 (r +1) (r ) + A 22 (r +1 \t ) (55) 

are generated by numerical integration. For (y +1) st iteration, these solutions are 
denoted by Zi H2 < m H1 ', • • • • rh Hn . The boundary conditions 

used in generating these solutions are 


r 0 ) = £, 
x i H Kt 0 ) = Q_, 



(56) 

i" n (r 0 ) = <l, = [OOO ... 1 ] r . 

Next, one particular solution at (y+1) denoted by xj , rk p , is generated by 
numerically integrating equation (52) from t 0 , to tf , using the boundary condi- 
tions Xi P (t 0 ) = Xo . Qi P (,t 0 ) — fl. Then, the complete solution of equation (52) can 
be obtained by using the principle of superposition and is of the form 


^U+D( r ) a- C - [ x 4 Hl {t ) + c 2 x i H Kt ) + .... + C n Xi Hn i.t ) + x i p (t ) (57) 

2h 0+1) (r ) = c lIU H Kt ) + C2 Tk" 2 (r ) + .... + c n ru Hn (t ) + ) (58) 

where the values of c x , c 2 c n which make 2k ( ‘ +1 0y ) = H/ are to be 

determined. To find these values of c x , c c„ , we let t = t f in equation 
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(58) and write it as 


ZV = 


21 i H Ktf) 21 i H Ktf) ru Hn {t f ) 


l-i 


Here, c_ A. 


ci c 2 


c r 

• • • • 


c_ + 21 i p (tf ). 
is unknown. Solving for c. yields 


(59) 


c = 


lu H Kt ) 21 (r, ) ru Hn (r, ) 




l-i 


a/ 


(60) 


It is important to note that the indicated matrix inversion in equation (60) 
has to exist in order to solve for c.. Substituting this solution of c_ into equations 
(57) and (58) gives the ( y + 1 ) st trajectory. This completes one iteration of the 
quasilinearization algorithm and this trajectory can be further utilized to begin 
another iteration, if required. Generally, the iterative scheme is terminated by 
comparing the y th and j + 1 st trajectories by calculating the norm shown in the 
following equation and comparing it with a preselected termination constatnt, p. 


Closed Loop Control 


h 

m 




<p 


(61) 


In order to obtain the closed loop control, the solution of the linearized equa- 
tion (47 ) can be written as 


X( */ ) = <t> ( tf ,t )x (r ) + f <t> ( t f ,r) [f (r) — / x ( r ) 1^ T (62) 

t 1 1 

where <t> is the state transition matrix of the system in equation (47). Rewriting 
equation (62) in terms of solutions of states and costates and replacing the integral 
terms by (t ) 
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*i (?/ ) 


011 f ) 012 (f/ » ? ) 

i (?) 


Hii C?) 

Hi (?/ ) 


021 (*/»*) 022 (*/ > ? ) 

Hi (?) 

+ 

Hi 2 C? ) 


From equations (42) and (63) 

m Cr 7 ) = = 0 2 i (r ; , ? )*j (?/ , r ) + 0 22 (r 7 , r )a (? ) + 


Thus, 


2L (? ) = 022 1 


“--02iSi Cr)-/> t2 (f) 

Oi 


= 022 1 


8 fi 

a & 


-£i2 (?) 


— 022 1 021 (? ) 


It is important to note here that 0 2 2 -1 always exists since it is 
minor of the state transition matrix. 

Substituting equation (65) into equation (38) 


(63) 
(?) . (64) 

(65) 
a principal 


Mi =£i (i. t ) . 


( 66 ) 
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APPENDIX A 

The nonlinear term N z is given in terms of the attitude of the orbiter, 0 _, and 
angular velocity vector transformation matrix, M , as 


N 2 = M" 1 


’A/" 

00 i 

'r 1 -^ 
00 2 

" . - — 1 r\M 


- M\ 


7 0 oi + A 22.1 


(A-l) 


where 


M = 


cos0 2 cos0 3 — cos0 2 sin0 3 sin0 2 
sin0 3 cos0 3 0 

0 0 1 


(A-2) 


o> = M t 0 


(A-3) 


and 


= [ o 0 0 1 
00i l I 


CA-4) 


T _ 1 

“ 00 2 COS0 2 


(-6>jSin0 jcos 2 © 3+to 2 sin0 2 sin 0 3 cos 0 3 ) (o^sin© 2 sin 0 3 cos 0 3 


— o> 2 sin0 2 sin 2 0 3 ) (t^cos© 2 cos0 3 — o> 2 cos0 2 sin0 3 )] CA-5) 


I (o> 2 cos0 2 ) ( OjCos© 2 ) 0 1. 

Q0 ? cos0 2 l 1 


(A-6) 


Since the transformation matrix, M , is a function of 0 2 and 0 3 , the time 


derivative of M can be expressed by the chain rule as 
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M 


4 = -2^-e, + 


002 Z 00; 


(A-7) 


From equation CA-1) 




(— sin0 2 cos0 3 )0 2 (sin0 2 sin0 3 )0 2 (cos0 2 )0 2 

0 0 0 

0 0 0 


(A-8) 




(— cos0 2 sin0 3 )0 3 (— cos0 2 cos0 3 )0 3 0 
(cos0 3 )0 3 (— sin0 3 )0 3 0 

0 0 0 


Substituting these equations (A-8) and (A-9) in (A-7) 


(A-9) 


M = 


(— sin0 2 cos0 3 )0 2 +(— cos 0 2 sin0 3 )0 3 (sin0 2 sin0 3 )0 2 +(— cos0 2 cos0 3 )0 3 (cos0 2 )0 2 
(cos0 3 )0 3 (— sin0 3 )0 3 0 

0 0 0 


(A-10) 


From equation (A-3), this can also be expressed as 


M = 


1 

COS0 2 


(— sin0 2 cos0 3 )(o> 1 cos0 2 sin0 3 
+6> 2 cos0 2 cos0 3 )+(— cos0 2 sin0 3 ) 
(— ojjsin© 2 cos0 3 +6i 2 sin0 2 sin0 3 
+co 3 cos0 2 ) 

0 

(cos 0 3 )(— ojjsin© 2 cos0 3 
+<o 2 sin0 2 sin0 3 +o> 3 cos0 2 ) 

0 

0 

0 


(sin0 2 sin0 3 )(o) 1 cos0 2 sin0 3 
+O) 2 COS0 2 cos0 3 )+(— cos0 2 cos0 3 ) 
(— o^sin© 2 cos0 3 +6> 2 sin0 2 sin0 3 
+O) 3 COS0 2 ) 

0 

(— sin0 3 )(— cOiSin© 2 cos0 3 
+6> 2 sin0 2 sin0 3 +co 3 cos0 2 ) 

0 

0 

0 
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cos0 2 (^icos0 2 sin0 3 +co 2 cos0 2 cos0 3 ) 
0 
0 
0 
0 
0 
0 
0 
0 
0 

Also, M~ l is given as 


(A-l 1 ) 




1 

COS0 2 


cos0 3 cos0 2 sin0 3 
— sin0 3 cos0 2 cos0 3 
0 0 


— sin0 2 cos0 3 
sin0 2 sin0 3 
cos0 2 


Thus, the nonlinear term N 2 can be rewritten as 


(A-12) 


N 2 = A' 3 (<* j 6J 


I 0 o)+A 22. 


Where the term A ' 3 is 


A ' 3 ((J,0J = M 1 


0 


802 

- M 

1“ a0 3] 



N 2 = A ' 3 (w,0J/ o co + A ' 3 (<±£JA zi 


= A 4 ((o,0J + A 5 (a0J)i. ( A- 1 3 ) 

where A 4 depends on the rigid-body slewing and is nonlinear in terms of and 0_. 
The second term relates the coupling between the rigid-body slewing and the flexi- 
ble modes. 
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APPENDIX B 

The transf ormation that relates the orientation angles 0. to Euler parameters j_ 
is a nonlinear transformation. This transformation is developed for body-three 
angles representation in this appendix and similar transformations can be derived 
for other three representations, namely space-three angles, space-two angles, and 
body-two angles. 

(a) For sin0 2 ^ 1 : 

If -JL <0 < JL then 

2 L 2 


0 j = sin 1 


2(^3 €i + €264) 


If ( cos0 1 cos0 2 ) ^0. then 


(B-l) 


0 1 = sin 1 


-2 U 2 e 3 - €i « 4 ) 

COS 

sin 1 

2 ( € 3 €i + € 2 €4 ) 



If ( cos0 1 cos0 2 ) <0, then 

0 ! = n — sin -1 

If C cos0 2 cos0 5) ^ 0, then 



-2 ( 

COS 

sin -1 1 


0 » = sin 1 


—2 ( €1 e 2 € 3 € 4 ) 


cos 


sin 1 2( e 3 €! + c 3 e 2 ) 


(B-2) 


CB-3) 


(B-4) 


If ( cos0 2 cos0 3 ) <0, then 
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2 C €i €2 €3^4 ) 

cos sin -1 2( €3*! + £3 £ 2 ) 

(b) For sin 0 2 = * 1 . 9 2 I s a constant. For sin 0 2 =l, 0 2 =-^-. However, if 

sin0 2 = — 1, then 0 2 = — j-. For this case, if ( sin0 isin0 2 sin0 3 + cos0 3 cos0 j ) ^0, 
then 

0 j = sin 1 2 ( e 2 €3 ■+■ €j €4 ) . (B-6) 

If ( sin 0 jSinO 2 sin 0 3 + cos 0 3 cos 0 3 ) < 0 , then 

0 i = it — sin -1 [2 C €2 £3 H- €3 €4 ) 1 . (B- 7 ) 


(B-5) 



For this entire case, 0 3 = 0 . 
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APPENDIX C 


Numerical Data: 

The analytics developed in reference [l] are utilized together with the basic 
SCOLE data [2] and the three dimensional linear vibration analysis [25] to generate 
the following numerical data. 


r = 


0.036 

-0.036 

-0.379 


Io 


1216640 -1.530307 175667.1 

-31.66433 7082976 -52474.84 

175690 -52503.9 7131493 
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0.45879£+2 

0.36305£— 1 

— 0.89042£— 1 

— 0.14067£0 

0.145 7£0 

0.36305£— 1 

0.62 11£ +2 

0.11263£0 

— 0.1471£0 - 

-0.5518£— 1 

-0.89042£— 1 

0.11263£0 

0.3273 7£ +2 

— 0.6392£ — 1 

— 0.14526£0 

— 0.1406 7£0 

— 0.1471£0 

— 0.6392£— 1 

0.254 7£ +3 

0.1908£0 

— 0.1457£0 

— 0.5518£— 1 

— 0.14526£0 

0.1908£0 

0.8103£+3 

0.1914£— 1 

0.19839£ — 1 

0.7925£— 2 

— 0.4278£ — 1 ■ 

-0.2570£ — 1 

0.84597£ — 1 

0.3935£— 2 

— 0.8369£— 1 

-0.761 15£-1 

— 0.12912£0 

— 0.6893£— 2 

— 0.7165£— 2 

— 0.2829£— 2 

0.1543£ — 1 

0.9222£— 2 

— 0.4269£— 1 

0.5969£— 2 

0.89767£— 1 

0.285 9 £—l 

0.461 1£-1 

0.4204£— 2 

0.41227£— 2 

0.1866£— 2 

-0.906 7£ -2 • 

-0.5 94 7£ -2 

0.1914£ — 1 

0.84597£— 1 

-0.68 93£ -2 

— 0.4269£ — 1 

0.4204£— 2 

0.19839£ — 1 

0.39 35£ -2 

— 0.7165£— 2 

0.5%9£— 2 

0.4127£— 2 

0.7925£— 2 

— 0.8369£— 1 

— 0.2829£— 2 

0.8976 7 £-l 

0.1866£— 2 

-0.4278£ — 1 

-0.761 15£— 1 

0.1543£— 1 

0.285 9 £-l 

-0.906 7£ -2 

-0.2570£ — 1 

— 0.12912£0 

0.9222£— 2 

0.461 1£-1 

-0.594 7£ -2 

0.23209£+5 

0.10383£— 1 

— 0.2089£— 2 

— 0.3955£— 2 

0.1227£— 2 

0.10383£— 1 

0.55561£+5 

— 0.37286£— 2 

— 0.3859£— 1 

0.2397£— 2 

-0.2089£— 2 

— 0.37286£— 2 

0.1342962£+8 

0.142 1£— 2 

-0.442 7£ -3 

— 0.3955£— 2 

— 0.3859£— 1 

0.1421£— 2 

0.2095672£+8 

— 0.9108£— 3 

0.1227£— 2 

0.23 97£ -2 

— 0.4427£— 3 

— 0.9108£— 3 

0.8662547£ + l 0 


Al 


— 0.2133821 E 0 
0.3808921£+3 
-0.1808478£+3 
0.1423380£+3 
-0.2416743 E +2 
— 0.6802273£0 
0.2784792£+2 
0.78428 18£ + 1 
-0.2694455£+2 
— 0.9225328£— 1 


-0.368705 7£ +3 
— 0.3030935£+2 
— 0.1318596£+3 
—0.113585 1£ + 1 
0.574383£+2 
0.3104929£2 
0.6651585£+2 
—0. 19 3009 7£ +2 
-0.5544 25 2£ +2 
0.1594045£+2 


•0.725390 1£-1 
-0.842765 8 £-l 
— 0.125799£0 
■0.236735 1£-1 
— 0.9150328£ — 1 
— 0.384306 2£ — l 
0.5 9607 5 £—l 
— 0.4363533£— 2 
—0.420062 3£ — 1 
— 0.1626004£— 1 
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The stiffness matrix K is calculated using equation (57) and the mode shape 
coefficients given in Table 1. This matrix is a diagonal matrix and is represented in 
terms of the diagonal elements as 


K = 


* u = 0.28202 17E0 
* 2,2 = 0.35 74692E0 
* 3i3 = 0.2412807E 1 
* 4,4 = 0.5285 116E1 
* 3 ,5 = 0.1588654E2 
* 6,6 = 0.8573860E2 
* 7 , 7 = 0.1146118E3 
* 8 ,8 = 0.56 86 10 1£ 3 
*9 9 = 0.6254598E3 
* 10,10 = 0.2114612E4 


The damping matrix B used for this analysis is a diagonal matrix and for 
damping ratio £ = 0.003, it is calculated to be 


B = 


5 U = 0.9685964E— 3 
ft 2 2 = 0.1088608E— 2 
ft 33 = 0.2834016E— 2 
* 44 = 0.4256808E— 2 
ft 5>5 = 0.7387177E— 2 
ft 6>6 = 0.1719014E— 1 
ft 7>7 = 0.1984237E— 1 
ft 8,8 = 0.442 1234E-1 
ft 99 = 0.4633434E— 1 
ft 10,10 = 0.8527647E— 1 
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Fy»are 1- Position Vectors in Inertial Frame 
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d(-L,t) 



Figure 2 ~ Vectors in Body-fixed Frame 



